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Abstract 

We consider an equation with drift and either critical or supercritical fractional diffusion. 
Under a regularity assumption for the vector field that is marginally stronger than what is re- 
quired for Holder continuity of the solutions, we prove that the solution becomes immediately 
differentiable with Holder continuous derivatives. Therefore, the solutions to the equation are 
classical. 

1 Introduction 

The purpose of this paper is to study the regularity assumptions needed on a vector field b for the 
solutions to an equation with drift and fractional diffusion 

ut + b-\7u+{-Ayu = f (1.1) 

to be differentiable (and therefore classical) . Wc do not assume the vector field b to be divergence 
free. 

Wc will concentrate in the so called supercritical case s G (0, 1/2]. The case s = 1 is usually 
referred to as the critical case, and s G (0, 1/2) as the supercritical regime because the drift term 
is of higher order than the diffusion term. In [17] it is shown that if 6 G C^~^'' in space for s < 1/2 
or b G L°° for s = 1/2, then u becomes Holder continuous for positive time. In this paper we 
prove that with the slightly better regularity assumption b G (^i-^s+a j^j. ^ g (0,2s), we 
obtain that the solution u becomes C^'" in space. This is large jump in regularity for a seemingly 
minimal extra regularity assumption in b. Note that b G C^~^* is the assumption that matches 
the scaling of the equation. A slightly better assumption like b G C^-'^^+°' allows us to use local 
perturbative techniques and thus obtain much better regularity results. 

For the case s > 1/2, a related result has recently been obtained in ]T5|, where a similar 
regularity for an eigenvalue problem is established. 

Our main result is the following. 

Theorem 1.1. Letb be a vector field in L°°{[~l,0],C^~^''+"{Bi)) and f G L°°{[~l,0],C^-^''+°'{Bi)) 
for some a G (0,2s), then any bounded solution u of (jl.ip in [—1,0] x Bi is C^'" in space. More- 
over, an estimate holds 

I'"IIl=°([-1/2,0],C1.°(Bi/2)) - ^ (ll"IU~([-l,0]xR") + II./IIl'=([-1,0],C1-2='+°(Bi))) 

where the constant C depends on s, n and ||&||pi-2s+Q only. 

If 6 G C^~^* and / G C^~^'* (case a = 0) we can prove that u is almost Lipschitz in space 
making a mild local smallness assumption. We state that as a second theorem. The proof is 
contained in the proof of Theorem 11.11 



1 



Theorem 1.2. There exists a S > such that the following result holds. Let b G [—1,0] x Bi be 
a vector field such that for some r > and any t G [—1,0] , 



sup 

\x-y\<r 



h{t,x)-b{t,y) 



< 5- 



(1.2) 



Assume also that f £ in space, then any bounded solution u of (jl.ip in [—1,0] x Bi is 

in space for all j3 < 1 . Moreover, an estimate holds 



where the constant C depends on r, s, n, /3 and ||6| [(71-23 only. 

The focus of this paper is on the regularity estimates for the solution u and not on the existence 
or uniqueness for a given vector field b and initial value u{x,0) = uq. It is our intention to use 
this result as a useful criteria to conclude that a solution is classical that could be applied to a 
variety of nonlinear equations with drift and fractional diffusion. Some applications are: 

• It implies the result in [5] about a regularity criteria for the supercritical surface quasi- 
geostrophic equation. 

• It complements [4] and |12| to prove that the solutions to the supercritical Burgers equation 
become classical for large time, instead of only Holder continuous. 

• Combined with |17j . it implies the result in j6l. Moreover, it proves a more general version 
of the result where u = L(— A)~*0 for any operator L of order zero (any combination of 
Riesz transforms or singular integrals), without assuming that divu = 0. 

• For any conservation law with fractional diffusion 



with F being a smooth nonlinear function, it proves that if 6* G 2s+a^ then actually 
9 G C^'" and it is a classical solution. 

• For the Hamilton Jacobi equation with fractional diffusion 

ut + H{Vu) + (-A)'u = 0, 

if H is smooth, it implies that if u G (71'1-2s+q some a > then, applying Theorem 
11.11 to the directional derivatives of u, m G C^'" for all a < 2s. In the case of s = 1/2, this 
provides an improvement of the regularity obtained in |18| . 

We also plan to use this result in future work for obtaining partial regularity results for non- 
linear problems for which some decay can be proved. For example, we expect that the methods 
in this paper would allow us to improve the results from [TB] and [3], although these two results 
have recently been improved with different methods in [8] and jl2j . 

We stress that there is no assumption on the divergence of b for either Theorem 1 1.1 1 or Theorem 
11.21 This allows us to apply the result to a larger class of equations. Moreover, there is no energy 
inequality for equation (jl.ip and variational methods do not seem to be suitable for a proof. 

The idea of the proof is to approximate locally the solution to the drift-diffusion equation by a 
solution to a drift-less problem via a change of variables (which is equivalent to approximating the 
solution locally by the solution to a problem with constant drift). The drift-less problem has C°° 
solutions, so we use it to show that there is a plane from which the solution separates slowly. A 



'^[[L°={[-l/2,0],Cf (B1/2)) - ([[w[[l°=([-1,0]xK") + [[/[[l°°([-1,01,C1-2=(Bi))) 



9t + dWF{9) + {-Ay 9 = 
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precise estimate on this separation leads to the C^'" estimate. This general idea goes back to the 
classical result of Cordes [7] and Nirenberg [M] for second order elliptic equations. For stationary 
integro-differential equation of order larger than one, it was developed in great generality in [5] . 
There is a difficulty that arises when trying to apply this method to integral equations of order 
one or less. The problem is that the afRne functions that are used in the approximations have a 
linear growth at infinity that already makes the tails of the integrals divergent when s < 1/2. In 
this paper this difficulty is overcome by rewriting the equation using the extension from [1] and 
complementing the affinc function in the original variable x with an extra term in the variable y 
with appropriate homogeneity (depending on s). 

The organization of the paper is as follows. In section [5] we explain the extended problem. 
The value 1 + 2s appears as the order of the second term in the expansion of solutions to the 
extended problem, and that is the main reason of the constraint a < 2s in Theorem 11.11 even for 
/ = 0. In section |3] we discuss the applicability of Theorem 1 1.1 1 to viscosity solutions. Essentially, 
the proof uses only the comparison principle with smooth sub and supersolutions, and thus it 
adapts perfectly to the case of viscosity solutions. In sections |3] and [SJ we develop the necessary 
regularity estimates for the drift-less problem. These estimates are not difficult and they take a 
large proportion of this article. But they need to be proved because the fractional heat equation 
is not classical in its extended version. The proofs in section |4] and [5] use essentially the same 
ideas as the usual heat equation and hence they are very standard. On a first read of this paper, 
it would be convenient to just skim through the results in these two sections. In section [51 the 
Holder estimates from are adapted to the extended problem. Sections [71 [SI and [HI are where 
the most important work of this paper is carried out specifically in order to prove Theorem ll.il 
Notation: 

• The variable x will always be used for a point in M" and y to be the extra coordinate in M+ 
for the extension to the upper half space. Sometimes the capital letter X will be used to 
denote a point in the upper half space M" x R"*". Thus X = {x, y). 

• By we denote the half ball centered at (0, 0) in the upper half space, i.e. = {(x, y) G 

X R+ : Ixp + y2 < i}. 

• By Bi we denote the unit ball either in the original space M" or in the boundary of the 
upper half space ]R" x {0}. 

• By (0_Bi)+ we denote the upper half part of the boundary of the unit ball dBi for which 
y>0. Thus {dBi)+ = dBi D {y > 0}. Also dB+ = U {dBi)+. 

2 The extended problem 

In [1] , the fractional Laplacian was characterized as a Dirichlct to Neumann problem for a degen- 
erate elliptic equation. More precisely, given a function / £ M" K, we can compute (—A)*/ by 
solving the following Dirichlet problem in the upper half space 

u(x,0) = f(x) inR" 
^ ^ ^ ^ ^ (2.1) 

divy'^Vu = in M" x R+ 

Then (—A)*/ = climy_>o Here Vu stands for the full gradient respect to the original 

variable x G i?" and the extended variable y G R"*". The constant a is equal to 1 — 2s. The most 
effective way to remember the relation between a and s is to use dimensional analysis. The order 
of the operator (—A)'' is 2s, whereas y°'dyU is of order 1 — a {—a because of y° and -1-1 because of 
dy), thus 2s = 1 — a. 
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Note that since in this paper we focus on the case s e (0, 1/2), the value of a stays in the range 
a e (0,1). 

We rewrite the equation (jl.ip using the characterization of the fractional Laplacian as a Dirich- 
let to Neumann operator given in [T]. 

ut{t, X, 0) + b{t, x) ■ - c hm y°-dyu{t, x, y) = f{t, x) for x £ M" and t G (0, oo) 

y-'O (2.2) 
div Vu = for j/ > 

The multiphcative constant c depends on the value of s (or, equivalently, on a). In terms of 
the smoothness of the solution of (|2.2p . the value of c is irrelevant, since it can be modified with 
a simple change of variables. So we will study the equation (|2.2[) with c = 1. 

The usefulness of this construction comes from the fact that the equation has become local. 
There is no integro-differential operator in (|2.2p . so the usual difficulty of keeping track of errors 
coming from the tails of the integrals that one has to deal with when studying regularity of 
integro-differential equations disappears in (|2.2p . 

In the next lemmas, we will explore the behavior close to y = of the solution to (|2.ip when 
the original function / is smooth. 

It is useful to remember some special harmonic functions for the extension (j2.1[) . The following 
list plays the role of the first few harmonic polynomials. 

• For any constant vector A e M", the function u{x, y) — A-x solves (|2.ip and limy^o y'^dyU = 
0. 

• The function u{x,y) = jir^y^"" solves (|2.ip and Ihny^oy'^dyU = 1. 

• For any constant vector A E M", the function u{x,y) = jh^y^ "" A ■ x solves (|2.ip and 
limy_j.o y'^dyU ~ A - x. 

• The function u{x,y) = a;^ — j-^j/^ solves (|2.ip and liuiy^Q y'^dyU ~ 0. 

Functions of the form A-x+ jzr^y^^"" play the role of the linear harmonic functions. Note that 
the next harmonic polynomial would be of the form A-x jz^y^'"^, which has degree 2 — a = l-|-2s. 
The following lemmas show that a first order approximation of smooth functions / close to the 
boundary in the extension problem (j9.1l) has an error of order 1 + 2s. This is the reason why in 
our main theorem, even if the right hand side is zero, we do not obtain a regularity estimate of 
the solution u in the space C^'^'* or better. 

Lemma 2.1. Let f : R" — > M &e a bounded function, f G C^'^{Bi), and u he its extension to the 
upper half space: 

u{x,Q) ^ f{x) onW X {0} 
div Vii = in R" x (0, +oo) 

Then, for x £ B1/2, u has the following expansion 

u{x, y) ^ fix) + yi-"g(x) + Oiy^). (2.3) 

More precisely, the error in the expansion is bounded by Cy^(| |/| l^ooj-jj,!-, + ||/||ci.i(Si)) for a 
constant C depending only on a and dimension. 
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Proof. Note that since / is bounded, u is bounded as well from the maximum principle. Therefore, 
the estimate of the expansion is trivially true for y > 1. We need to check the case of y smaller 
than 1. 

We compute u in terms of / explicitly using the Poisson kernel from [T]. 



y 



u{x,y)^Cs j fix + z)- ^ 

(|z|^ + y^j 2 



d^ 



We consider g{x) = — c(— A)'*/ and we have 

fix) =Cs [ fix) 



(|z|2 +?/2)- 



yi-»<7(x) = ifix + z)~ fix))- 



,1-a 



|7i + l-a 



Therefore, we estimate the error in 

\uix,y)-fix) -yi-"g(a:)| = 



ifix + z)^ fix))y 



l-o 



i\zV + v^)- 



|n+l — a 



dz 



Since the kernel is even, the first order part of / in Bi /2 intergrates to zero and we can estimate 
the integral with the C^'^ norm of /. 



< C / (|z| xSi/2 +XE"\Si/2)y 



+ y"^) 2 \z\ 



n+l — a 



dz 



In order to estimate the value of this integral, wc split the domain of integration in two subdomains: 
For |z| > y/2, we can estimate the numerator in the kernel as 



i\z\' + y')^ <C\z\ 
For z G R" \ By/2 , we estimate the integral 

C I UIV"° _ " , dz = Cy^ 



V ■ 



(|z|2 +2/2)- 



dz 



Cy' 



•\B„/2 (|z/2/|2 + l)- 



— dz = Cy^ 



Finally, for |z| < y/2, we have 



|"+l-_(|^|2+y2):^ 



Cy 



n+1 — a 
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Thus, we estimate 



C I \zVy 



2„,l-a 



l"+l-'^-(|z|2 + y2)i±i^ 



|n+l— a 



Which finishes the estimate of the error in (|2.3p . □ 

The foUowing lemma is a locahzed version of Lemma 12.11 
Lemma 2.2. Let u solve 

w(x,0) = /(a;) on 
div y^Vu = m 

Assume that / G C^'^(i?2') mc? u G L^{B^ ,y°-). Then u has the foUowing expansion in B^j^. 

u{x, v) = f[x) + v^-^g{x) + 0(2/2) (2.4) 

for a function g G C^^"^ and the error is controlled by Gy'^{\\f\\c'^.i + ||M||L2(-j,a-)), where C depends 
on n and a only. 

Proof. Let v be the extension to the upper half space of the function / extended to R" as zero 
outside B^. 



f{x) for |a;| < 1 
for \x\ > 1 



v{x,0) 

div y'^ Vw ^ in M" x (0, +oo) 



We apply Lemma I^TT] and obtain an expansion like (|2.4p for v. 

We are left to show that the remainder w — u ~ v has an expansion of the same kind. From 
the maximum principle, \v\ is bounded by and therefore ||i'|li2(^+ j^a) < 011/111,=° and 

We can consider the odd extension of w to the full ball Bi by w{x, —y) = —w{x,y) so that w 
solves the degenerate elliptic equation across {y = 0}: 

div|y|°Vw = in Si. 

We apply the Harnack inequality from jllj to w. Then w is Holder continuous in the interior of 
Bi. Moreover, from the Cacciopoli inequality 

f y"|Vwp <C f y''\w\'^. (2.5) 

We will apply the boundary Harnack principle from jlOj to w in Bf. In order to apply that 
theorem, we split w into its positive and negative parts 



= on 9J5j^^ 



div y^VujP = in B^ 



/4 



= w on dBy^ 



divy"Vw" = in S, 



3/4 
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Thus, w — wP — w", wP — w" = on -B3/4 and both and w" are nonnegative. Since 
div y°V(?/^~") = in the upper half space {y > 0} we can apply the Boundary Harnack theorem 
from [TUj to obtain that both /y^~'^ and arc Holder continuous up the the boundary 

in S^g- Therefore, there exist a Holder continuous function g such that 



w{x, v) = .9(2;, y)y ° in B, 



5/8- 



Recalling that u = u + w, so far we have proved that u has an expansion of the form u{x, y) = 
f{x) + y^^°'g{x) + 0(2/^^"+") for some a > (the Holder exponent in g). 

We can repeat the same argument for any derivative (or incremental quotient) of w with 
respect to x. From Cacciopoli's inequality, all those derivatives are in L^{y"-) and also solve the 
same equation. Therefore we obtain 



Wx,ix,y) ^ gx,{x,y)y^ " in S+^g. 

with gxi Holder continuous. Repeating this argument, we obtain that g is C°° in the x variable. 
Now consider 7/2(2;, y) ~ y'^dyW. Then we have 

y-^'ul = / y^\dyw\' <C [ y^\w\' < C{\\f\\^^ + 

1/2(2;, 0) = 5(x,0) 

Moreover, as it is pointed out in [1], M2 satisfies the conjugate equation 

div?;"°Vu2 = in S+g 

So we can start over out argument with U2 instead of u and —a instead of a to obtain a Holder 
continuous function h such that 



U2{x,y) = g{x,0) + y +°7i(2;,y) in BJ^^. 

Therefore 

ry ry 
w{x, y)= dyw{x, z)dz= / dyy-^gix, 0) + yh{x, y) dz = y^^°-g{x, 0) + 0(y2) 

Jo JQ 

Which finishes the proof. □ 
Remark 2.3. Note that the function g{x) corresponds to the fractional Laplacian of /. 



0(2;) = hm — = hm y OyUix.y) 



Corollary 2.4. Let f : B'f M. be a C^'^ function and u : Bf R satisfy the extension PDE: 

u{x,Q)^f on Bl 
div(?/°Vu) = tn B+ 

where a > 0. Then there exists a number D and C > such that 

\u{x, y) - /(O) - V/(0) • X - Dy^-'^\ < C{x^ + z/ + {xly^-") (2.6) 

where C and an upper bound for \D\ depend on \\f\\c^{B°) '^'^d 11^11^2(5+ yay 
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Proof. From Lcnima l2.2[ u has an expansion of the form 

Since / G C^ K, fix) = /(O) + x ■ V/(0) + 0{\x\^). 

The function g is C^+°. Since a > 0, then g is Lipschitz. If we choose D ~ g{Q), we have 

g{x)y'--^g{0)y'-'^ + O{\x\y'-'^). 

Adding the three terms together, we finish the proof. □ 



3 Weak solutions 
3.1 Viscosity solutions 

Since we do not make any assumption on the divergence of b, integration by parts is comphcated 
and the distributional sense is not suitable for defining weak solutions of p.ip . Since we are 
assuming that b is continuous, we can use the concept of viscosity solutions developed originally 
by Crandall and Lions for Hamilton- Jacobi equations. The most straight forward way to define it 
for the integral equation p.ip is the following. 

We say that a continuous function u is a viscosity solution of (|l.ip if for every point {xq, ta) € 
(-1,0] X and every bounded function : [-l,to] x I^" such that 

1. (ys is smooth around (to,a;o)i 

2. (f>{to,XQ) = M(to,Xo), 

3. (j)>u (or respectively (f> < u) in [— l,to] x , 
then 

(j)tito,xo) +b(to,xo) ■ V(j){to,xo) + {-Ay(l)(to,xo) < f(to,xo) ( or > f{to,xo) respectively.) 

This is the definition used in [18]. in the same way as distributional solutions of PDEs are 
defined based on the property of integration by parts against smooth functions, viscosity solutions 
are defined based on the property of the comparison principle with smooth sub and super solutions. 
This concept is more apparent in the following definition of viscosity solution for the extended 
problem (j2.2p . Indeed, the definition says that a continuous function w is a viscosity solution if it 
satisfies the appropriate comparison conditions with respect to smooth sub and super solutions. 
We give the precise definition below. 

We say that a continuous function u : [—1,0] x K solves p.2|) in the viscosity sense if 

for any t S [—1,0], u{t,x,y) solves the Dirichlet problem (|2.ip (classically) in {t} x B^ and for 
every function (p : [to — e, to] x B^{xo) — R such that 

1. ip is smooth in [to — £,to] x ^^(a^o), 

2. ip satisfies divy^Vt/) = in B^{xq) for all t <E [to — e,to], 

3. (fit + b ■ V0 — limy^o y'^dycf) > / (or < / respectively) in (to — e, to] x B^{xo) 

4. (f> > u {or < u respectively) on {to — e} x Bf{xo) U [to — e,to] x {dB^{xn))^ , 
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then (j) > u or (j) < u respectively, in the whole domain [to — e, to] x Bf{xo). 

It is not hard to prove that the definitions of viscosity solution for ([OJ implies the correspond- 
ing definition for (|2.2p after the extension (|2.ip . 

Viscosity solutions is the most appropriate type of weak solutions for the results in this paper. 
The actual equation is used directly only in the proof of Lemma 17. 1[ and the way it is used is to 
compare the solution with a classical subsolution. 

Note that for any other type of weak solution u (for example entropy solutions), the same 
function would solve the equation in the viscosity sense as long as the comparison principle with 
classical sub and super-solutions holds. In practical cases, a drift-fractional diffusion equation may 
be obtained as the linearization of an equation for which the comparison principle can be shown 
to apply, and thus the results of this paper would apply as well. 



3.2 The vanishing viscosity method 

One way to avoid the concept of weak solutions completely is by adding a vanishing viscosity term 
and proving uniform estimates. This method is usually called the vanishing viscosity method. In 
this context, it is not related to the concept of viscosity solutions at all despite the unfortunate 
similarity of the names. 

The idea consists in adding an artificial viscosity term to the equation 

ut + b- Vu + {-Ayu - eAu = / in [-1,0] x M". 

With the added artificial viscosity eAu, the equation has a smooth classical solution. 

If we prove that the estimates in this paper hold uniformly as e — 0, then the existence of 
classical solutions for (jl.ip follows by passage to the limit. In terms of the extended equation, the 
same idea applies to 

ut + b-\7u- {\im y°^dyu) - eAu = f in (-1,0] xS^*. 

divy''\7u = Q in [-1,0] xB+. 

The estimates in this paper are indeed uniform as e — > 0. The explicit form of the equation 
(|2.2p is used only for Proposition 16.11 and Lemma 17.11 The first is a localized version of the result 
in |17| . which adapts to the vanishing viscosity method (as pointed out in that paper). The second 
also adapts to the vanishing viscosity method as it is analyzed in Remark 17.31 

The vanishing viscosity method is a very common approach for nonlinear conservation laws. 
It is used for example in [3], [13] and [9] for problems with drift and fractional diffusion. 

In this paper we prefer to focus on the estimates but we will point out precisely how the proofs 
adapt to either continuous viscosity solutions or to the vanishing viscosity approximation. It turns 
out that the precise form of the equation is used only in very specific steps in the paper, so it is 
only necessary to analyze the type of solution at those points (See remarks 17.21 and 17. 3p 



4 Interior estimates for the fractional heat equation 

In order to carry out the perturbative arguments to prove the main theorem, we need to have 
good local regularity estimates for the drift-less equation. 

ut{t, X, 0) - lim y°-dyu{t, a:, y) = on (-1, 0] x i?" 

divy'^Vw^O in[-l,0]xBj^ 
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For any time t e (0,1), the solution of (|4.ip is going to be C°° in in the interior of B^. 
Moreover, u and all its partial derivatives of any order respect to x are Lipschitz in time. The 
smoothness with respect to y can then be understood by Lemma 12.21 

Proposition 4.1. Let u be a bounded continuous solution of the equation (|4.ip . then u and ut 
are C°° with respect to the x variable in the interior of the domain. Moreover, for any multiindex 
a, the following estimate holds 

II^""IIl~([-1/2,0]xB+2) - II"IIl«=([-1,0]xB + ) ' 
||a"9tu||L~([-l/2,0]xS0^2) - II'"IIl°-([-1,0]xS + ) ■ 

Here the constant Ca depends on a, s and n only. 

Proof. Note that if we show an estimate for u in C^{By^), then the estimate in -8^2 follows 
by the fact that, for each time, m is a bounded solution of the Dirichlct problem for the elliptic 
equation of the extension. The same argument cannot be made for ut since a priori we cannot 
conclude that Ut is bounded in the extension where y > 0. 

By a classical translation and scaling argument, it is enough to show the estimate at the point 
(0,0,0). 

Let h{t, x) be the fractional heat kernel, 

ht + i-Ayh = in (0, +00) X M" 
h{x, 0) = Sq 

There is no known closed form for h{t,x) in real variables. However, in Fourier variables it 
is simply h{t,£^) ~ e^*'^' . Let us extend h to the upper half space h{t,x,y). Using the scale 
invariance of the equation, we have that for some function H : R" x — > , 

h{t,x,y) =t"^H (t^^x,t^^y^ (4.2) 

The function H is the extension to the upper half space to the function whose Fourier transform 
is e~'^' . It is bounded, smooth in x, and H{x,0) ~ (1 + |a;|)~"~^*. From the Poisson formula in 
[T] we also obtain \D^H{x, y)\ < C(l + yY^il + |a;| + y)-"-2.'s-fc 'pj^^, function h is bounded and 
C°° in X and t except around the point (0, 0, 0). 

This function h is the heat kernel of fractional diffusion in the whole space. We will multiply it 
by a smooth cutoff function in order to find an expression of w(.t, y, t) in terms of localized integral 
quantities. Let 77 be a smooth radially symmetric function in x and y, supported in By^ and 

identically one in B'^^^. Note that since i] is radially symmetric and smooth, dy-qix^y) = 0{\y\) as 
y ^ 0. In particular limj^_>o y°'dyr]{x, y) = 0. 

At any point x G '^^^ recover u(t,x) by the formula 

u{x,t) = lim / u{t, z)h{t ~ T, X — z)ri{x — z) dz. (4-3) 
•^^0+ J go 

This formula allows us to compute estimates for solutions u of (j4.ip in a bounded domain. 
Let h(t,x,y) = h{t,x,y)r]{x,y). Note that since dyr]{x,y) = 0{y) near {y = 0}, h satisfies 

ht - lim y"^dyh = 0. 
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and also divy^V/i is a bounded function supported in \ ^3/4 which is C°° in x. Let us 

analyze its behavior as t — ;> 0. 

div y°-Vh{t, x, y) = div (y" V [hrj)) 

= h div (y" V77) + V/i • V77 

= f-^Hif-i^x, t-i^y) div (y^V?/) + 2yT^^^V H{r^ x, r^y) ■ Vry 

The last estimate uses that H{x,y) « {\x\ +?/)^"^^'* and VH{x,y) w (|a;| + y)^"^^''^"^ for large 
\x\ + y. 

Using (|4.3p . we compute 

w(0,x,0)— / u{—l,z)h{l,x~z)dz^ / 9t / t, a; — ^) d^ dt 

Jb? J-1 Jb'I 

/ Wt(t, z) h{—t, X — z) — u{t, z) ht{—t, X — z) dz dt, 
-1 Jb° 


-{-l\yu{t,z) h{-t,x~ z) 



1 J B 



-1 



+ u{t, z) {~AYh{-t, X — z) rj{x — z) dz dt 
u{t^z^y) div {y'^V {h{~t,x ~ z,y)^^ dz dy dt 



Now we have u(Q,x, 0) written as a sum of two convolutions with functions which are C°° in 
X and compactly supported. Therefore u is C°° in x and the first estimate of the Lemma follows. 
From the equation (|4.ip . we have that 

Ut — lim y°'dyU. 

Since u is C°° in x, the right hand side is well defined and C°° in x. From this we deduce the 
second estimate of the Lemma. 

□ 

Remark 4.2. It is impossible to find a better regularity estimate in time in terms of the L°° 
norm of u only. The problem is that the estimate for dtu on cannot be extended to B^^^. 
One can show that if m = on (0, 1/2) x {dBi)~^ and m = 1 on (1/2, 1) x {dBi)~^ , then indeed 
the solutions to (j4.ip will be only Lipschitz in time on {1/2} x By 2- However, just assuming 
that u is continuous on (0, 1) x (dBi)^ would imply that u is in time from the expression of 
Ut ~ lim-y_j.o y°'dyU as a convolution. 



5 Boundary behavior for the fractional heat equation 

In this section we study the Holder continuity of the solution to the fractional heat equation (|4.ip 
on the parabolic boundary. All the estimates in this section are proved by comparing the actual 
solution with explicit barriers. The results are very natural. They simply say that if the Dirichlet 
boundary data for the fractional heat equation is Holder continuous, then the solution is Holder 



11 



continuous on the boundary. Wc tried to make the computations of the barriers as simple as 
possible, but some tedium is unavoidable. On a first reading of this paper, it would be advisable 
to skim through the results in this section. 

The estimates in this section follow from the local properties of the fractional heat equation 
(|4.ip . The equation is properly a parabolic evolution equation on [—1, 0] x 5°. For each fixed value 
of t, u solves the elliptic equation divy^Vu = on {t\ x . We can expect to find regularity 
estimates on [—1,0] x B^ and also on {t} x for each t. However, in this section we do not 
intend to analyze a modulus of continuity with respect to t for x, y) for positive values of y, 
since intuitively those quantities arc not connected by the equation. 

We start by constructing a barrier that will be useful to understand the continuity of u close 
to dBi for each fixed t. 

Lemma 5.1 (Barrier for ((9i?i)+). Let a G (0, 1) and — (xq, yo) be any point on dB^\{y — 0}. 
There exists a constant C > such that the function 

B{X) C(l - iXp)" + \X - Xol" 

is a supersolution of 

diYy°-VB<Q inB+. 
Proof. One way to simplify the computations is to note that 

divy'^VB ^ J/" (^AB + ^dyB 

Let us use this formula to estimate the equation in the two terms of the definition of B. We 
compute 

dWy'^V ((1 ~ \X\^r) = [a(1 - \X\'r + -dy{l - \X\^r 

= y" {-a2{n + a){l - iXp)""! + a{a - 1)4|X|2(1 - \X\^Y'-^) . 
= y'^ail - (_2(n + a)(l - \Xf) + \\X\\a - 1)) 

This is a negative value since a < 1. Moreover, it becomes very negative as \X\ gets close to 1. 
Now we estimate the second term. 

divy'^V(|X-Xor) = ay" (^n + a - 1 + a^— iX-Xol"-' 
Since yo is positive and a > 0, 

< (n + a - 1 + a) |X - XqI""^ 
For all X ^B'l,{\- |Xp)"-2 ^\^^ _ Xo|"~^ so we can choose a constant C > such that 

divy'^V(C(l-|X|r + |X-Xor) <0 
for all X e Bi+. □ 
Lemma 5.2 (Continuity on (9i?i) + ). Let u : B^ K. be the solution to 

div Vw ==0 in B^ 
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Assume that for some (x'OiJ/o) G {dBi)^ and every {x,y) G dB^ we have 
\u{xo, vo) - u{x, y)\ < Co((.T - xof + (y - yof)" 

Then for some C > 0, 

\u{xo,yo) - u{x, y)\ < CCo{{x - xof + (y - yofT 
for all {x,y) G . Here C depends on s, n and a. 

Proof. From the assumption, the function U{x,y) — w(a;o,yo) + C'oB(x,y) is larger or equal than 
u on dB^ , where B{x, y) = B{X) is the supersolution constructed in Lemma [01 Therefore, from 
comparison principle u <U in the whole half ball B^ . Therefore 

u{x,y) - u{xo,yo) < CoB{x,y) < CoC{\x - xo\ + \y - yo\)°'- 

□ 

The next barrier is useful to analyze the Holder continuity close to B^ of solutions to div y"'Vu = 
in B^ . This barrier is not used in the rest of this section. 

Lemma 5.3 (Barrier for B^). Let a G (0, 1 — a) and xq G i?". There is a C > I such that the 
function 

B{x,y) = \X -{xo, 0)1" + Cy" 
is a supersolution in the upper half space 

divy"VB<0 /or yG (0,1). 

The proof is a direct computation similar to the one of Lemma [5.11 so we omit it. 
We now construct an auxiliary function that will be useful to construct barriers. Let B be the 
only bounded solution of the following problem. 

B{x,y)^0 on{dBi)+ 
- \imy''B{x,y) ^ 1 in B*^ 

divy''VB = in B+ 

The explicit formula for B can be computed using the intuition from [1] that the equation in 
the upper half plane can be understood as the Laplace equation in fractional dimension. 

B{X)= [ $(z-X)- |X|-"-"+i$(z-X/|Xndz (5.1) 

where $(X) = l/|X|"^^+° for any X G M"+^ is the fundamental solution at the origin. It is 
elementary to check that this is the correct formula for B. 

From the formula (jS.ip we can see a couple of elementary properties of B: 



1. There exists c > such that B{x, 0) > c(l — \x\y . 

2. The maximum of B is achieved at x = y = 0. Moreover there exists c > such that 
B{0, 0) - B{x, y) > c(|a;|2 + y^""). 
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Lemma 5.4 (Continuity on dB^ x ( — 1,0)). Let u be a solution to the localized fractional heat 
equation (|4.ip . Assume that for some xq G dB^ and tg G [—1,0] we know 

\u{t, X, y) - u{to, xo, 0)1 < Co {\x - xo\^ + y'"" + (to - t))" 

for all {t,x,y) G ({ — 1} x B^) U ([— l,to] x {dBi)'^). Then there is a constant C > such that 

\uit,x,y)-u{to,xo,0)\ < CCo(|x-xoP+?/-'^ + (to-<))'"^' for all (<,x,y) G [-l,to] x B+ . 
Here C depends on a, n and a but not on to. 

Proof. Let r = {\x — a;op/4 + y^^"" j (1 — a) + (to — 0). Let U be the following barrier function 

U {t, x, y) = (xo -x)-XQ + -^y^'" + {to - t) + 2B{x, y) 

1 — a 

where B{x,y) is the function in ()5.1|) . 

We see that r <U < Cr'/'^ in [-l,to] x B^ . 

For a given {ti,xi,yi) G [— l,<o] x B^ , let us choose p = 4rJ^^ (where ri = \xi ~ a;op/4 + 
yi^^/ll - a) + (to-ti)) so that on [-l,io] x {dBi)+ U {-1} x B+ we have 

Cop"-iC/ + Cop" > u(t, y) - u{to, xo, 0) 

Since both the left hand side and right hand side solve the localized fractional heat equation, then 
also for all {t,x,y) G [— l,to] x B^ 

u{t, X, y) - u{to, xo, 0) < Cop'^-^U + Cop" 

In particular 

u{ti,xi,yi) - u{to, Xo, 0) << CCorl°'^^. 
The inequality from the other side follows by applying the same barrier bound to u(to, xo, 0)—u. □ 

Lemma 5.5 (Continuity on i = 0). Let u be a solution to the localized fractional heat equation 
(|4.1I) . Assume that for some xo G B^ we know 

\uit,x,y)-u{-l,xo,0)\ < Co(|a:;-xoP+yi-" + t+l)" (5.2) 

for all (t,x,y) G { — 1} x B^ U [—1,0] x dB^ . Then there is a constant C > such that 

\u{t, X, y) - u(-l, Xo, 0)1 < CCo ([a; - xqI^ + y^"" + t + l)""^^ for all [t, x, y) G [0, T] x B+. 
where C depends on a, n and a. 
Proof. Wc will use the following barrier function 

U{t,x,y) = \x~xo\^ + {-y^ + 2y^-'') + 2n{t + 1) 

(1 + a) ' 

which is a solution to (|4.ip . Let r = [a; — xol^+y^^^ + t + l. Note that there are constants c and 
C such that cr <U < Cr. 
From ((O)) . 

|M(t,a;,y) -u(-l,a;o,0)| < Cor" for all {t,x,y) G {-1} x U [-1,0] x {dBi)+ 
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Therefore, for any ri > 0, on the paraboUc boundary {t,x,y) G { — 1} x U [—1,0] x (dBi)'^ we 
have 

u{t,x,y)-u{-l,xo,0) < Cocr'^-'U + Cor'^ 

For any point (ti,Xi,yi) and corresponding value of ri, we apply comparison principle to obtain 
that 

u{t,xi,yi) - u(-l,a;o,0) < CoCr" . 
Using the same argument to bound u(— 1, xq, 0) — u we finish the proof. □ 

Lemma 5.6 (Main resuh of this section). Let u be a solution to the localized fractional heat 
equation ()4.1[) . Assume that ||w||^oc([_x g]xB+) 1^ C'o and the following regularity for the boundary 
data: 

Il"llc°({-i}xi3+) < (5.3) 
l'"llc°([-i,o]xas?) < C'o (5.4) 
ll"(i,-,-)llc"((OSO+) < ^0 for all te [-1,0] (5.5) 

There is a constant C such that if {t,x,y) is any point in [—1,0] x and {t* ,x* ,y*) is either 
its closest point on {t} x (9i?i)+ if y > or its closest point on ({ — 1} x Bf) U ([—1,0] x dB^) if 
y = 0, then 

\u{t,x,0) - u{t*,x*,y*)\ < CCa{\x - x*\ + \y - y*\ + \t - t*\f (5.6) 
for some a > depending on s, n and a only. 

Proof. We wiU start the proof by analyzing the case y = 0. In this case y* — y = 0. 

The point {t*,x*,0) belongs either to the bottom of the parabolic boundary: { — 1} x Bi or 
the lateral boundary [—1,0] x dB^. We will analyze each of the two cases. 

li {t*,x*,0) e {-1} X Bl, then from dO]), (EH) and ([53]), the hypothesis of Lemma [13] are 
satisfied (with a/2 instead of a). So we apply Lemma 15.51 and prove (|5.6p with a = s'^a/2 

If {t*,x*,0) e [-1,0] X dB^, then from dOJ, and ([53]), the hypothesis of Lemma [O] are 
satisfied (with a/2 instead of a). So we apply Lemma [5^ and prove ()5.6|) with a = s'^a/2 

So we have finished the proof in the case y = 0. What we proved so far implies that for every 
ito,xo,yo) G [-1,0] X {dBi)+ and {to,x,y) £ [-1,0] x dB+ (recaU dB+ = {dBi)+\JBl), we have 

\u{t\x, y) - u{t\x*,y*)\ < CCo{\x -x*\ + \y- y*\f (5.7) 

Let us analyze the case y > now. In this case {x* , y*) G dBi. On the other hand, from (15.51) 
and (|5.7p . we can apply Lemma 15.21 with xq — x* and y^ = y* and finish the proof. □ 

Remark 5.7. Note that the barriers in this section could be used (together with the comparison 
principle) to show the existence of the solution of the Dirichlet problem for (|4.ip by Perron's 
method. 



6 Holder estimates for the equation with drift 

Our assumption h G (7^^^''+" means that the vector field b is slightly better than C^^^'^. From 
[T7], just from 6 being in the class C"'^"^", we can obtain that the solution u is Holder continuous. 
In this section, we adapt the result of [17] to the extended problem which allows us to localize the 
estimate. 
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Proposition 6.1. Let u be a bounded function in [—1,0] x solving 
Ut + b-Vu- lim y°-dyu{t, x, y) = f on (-1, 0] x 

(6.1) 

div y^Vu = in [-1, 0] x 

where b is a vector field in [—1, 0] xBi which is C^"^" in x and f is a scalar function in L°°([— 1, 0] x 
Then u is Holder continuous in [—1/2, 0] x B^^^ '^^^ also in {t} x B^^^ for every t E [—1/2, 0] 

l'"llc°([-l/2,0]xS0^2) + II"IIl°°([-1/2.0].C°(S+2)) " "-^l I"' Il°°((-1,0] xB+) 

for some constants a > and C depending on s, ||6||pi_2s and the dimension n. 
Proof. Consider the function w : [— 1, 0] x R" x K+ defined by 



cc, 0) 

div y'^Vv = for y > 



u{t, X, 0) if |a;| < 1 
otherwise 



Note that the function u ~ v satisfies 

{v-u){t,x,0) =0 inB?, 
divy"V(w - u) = in B+. 

Moreover, ||u — it||^oo < 2 Thus, for every fixed vahie of t wc can apply Lcinma l2.2l to 

obtain that 

g{t,x) ^ \im y°-dy{v - u){t,x,y) 
is differentiable in x (with an estimate depending only on ||?i||^oe. In particular, g is bounded on 

[-1,0] Xi?30/4. 

Now, let us analyze the equation for v. We see that Vt = ut and Vt; = Vu on [—1,0] x B^. 
Therefore v satisfies the equation. 

Vy + 6 • Vw - lim y^dyvit, x,y) = f + g on (-1, 0] x B" 

div y^Vv = in [-1, 0] x 5+ 

Since the right hand side g is bounded in [—1,0] x -63^4, wc can apply the main theorem from |17) 
to obtain that v is Holder continuous in [—1/2, 0] x B^^. That means that u is Holder continuous 
on [—1/2,0] x i?2/3- From this, one concludes that u is Holder continuous on {t} x B^^^ for each 
t G [—1/2, 0] using the barrier from Lemma and that u{t, ., .) solves the equation divy^Vu ~ 
in B+. □ 



7 The approximation lemma 

The main lemma in this section says that if the vector field b and the right hand side / are very 
small, then u differs very little from the solution to the fractional heat equation with the same 
boundary condition. This lemma is important in order to carry out our proof by perturbation. 
Indeed, the main idea of the proof of Theorem 1 1 . 1 1 will be to approximate the solution u with the 
smooth solution of the drift-less problem with zero right hand side. 
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Lemma 7.1. Let u be a continuous Junction that solves 

ut + b ■ Vu - lim y°-dyu{t, x,y) = f on (—1, 0] x i?° 
div Vu = m [-1, 0] x B+ 

Assume that ||u||^oo < 1, l|w||(7c.(-[_i qjxbO) ^ C* and IKtllioo^.i qj ^0(5+)) < C for some a > 0. 
Let V be a solution of the driftless problem 

vt{t,x,0) - lim y''dyv{t,x,y) = on (-1,0] x B'^ 

divy^Vv^O ot[-1,0]xB+ (7-2) 
v{t,x,y) = u{t,x,y) on the parabolic boundary: ([—1,0] x (dBi)'^) U ({ — 1} x B^). 

Then for every e > there exists a S > so that ||w — u||^oc([_i o]xs+) ^ ^ every time that 
II^IIl~([-i,o]xbO) < S and 1 1/| |l'-([-i,o] xBf ) < ^- 

Proof. From Lemma [5.61 the for any {x,y,t) e [—1,0] x B^ and {x*,y*,t*) as in Lemma [5.61 we 
have 

Hx,y,t) - vix*,y\t*)\ < Ci\x - x*\ + \y ^ y*\ + \t - t*\f 
On the other hand, since ||u||p<:.Q„]^ o]xb°) — ^ II"IIl«=([-i 0] c°(b+)) — ^^^^ ^lave 

\uix,y,t) - uix*,y*,t*)\ < Ci\x - | + |y - | + |t - t*\f 

But v{x*,y*,t*) = u{x*,y*,t*). Therefore, if {x,y,t) befongs to either [-1,0] x (S+ \ B+_p) or 
[-1, -1 + (5] X Bl, we have 

\v{x,y,t) - u{x,y,t)\ < Cp"". 

Moreover, on each time shce t G [— 1, — 1 + (5], we use maximum principle for the degenerate elliptic 
equation 

divy''V(u - u) = on {t} x B^ . 

Then |u — u| < Cp" in [— 1, — 1 + (5] x B^ . Choosing p small enough, we can assure that |u — 1;| < e/2 
in([-l,0]xS+)\([-l + p,0]xB+_^). 

From Proposition 14.11 Vi; is bounded in [—1 + p, 0] x i??_p with a bound depending only on 
dimension, a, p and jji'lj^^oo = 1. Therefore 

vt{t,x,0) + b-\/v - lim y''dyv{t,x,y) < C6 in [-1 + p,0] x S?„„ 

We can conclude using maximum principle that v = v — e /2 — Cdt < u. Indeed, we know that 
-D < M on the boundary and complement of [—1 + p, 0] x Bi^p, whereas inside we have 

vt{t,x,0) + 6 • V5 - lim y'^dyv(t,x,y) < -S. 

y^O 

Similarly, we prove that v + e/2 + CSt > u. Wc finish the proof by choosing S small so that 
CS < e/2. □ 



Remark 7.2. Note that the fact that u solves the equation p.2p is used only to apply the 
comparison principle with the classical subsolution v, thus u only needs to satisfy p.2p in the 
viscosity sense. 
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Remark 7.3. The proof of Lemma [7.11 adapts with little work to the vanishing viscosity frame- 
work. Indeed, if the equation for u is 

Ut + h-S/u- (lim -eoAu = / on (-1,0] x 5° 

(7.3) 

divy^Vw^O in[-l,0]xBJ^ 

then under the same hypothesis plus £o < <5, the solution to (|7.2p approximates (|7.3p . The only 
modification in the proof is that when we analyze the equation for v we also have the term £o/\v, 
which is small since by Proposition 231 w is in x. 



8 The improvement of flatness lemma 

Our main regularity result is proved iteratively by showing that the values of the solution separate 
less and less from a plane in smaller scales. In each iteration we use the essential improvement of 
flatness lemma which is stated and proved below. 

Lemma 8.1. Let a he any positive number less than 2s, b G C^^^", and u be a continuous function 
that solves 

ut + b- Vu — lim y°'dyu{t, x,y) — f on (—1, 0] x 5° 

(8.1) 

divy'^Vw^O m[-l,0]xB+ 

Assume that ||w||^oo < 1 and ||6||^oo < S and < S. Then, if S is small enough, there is a 

radius r > 0, a vector A e M" and a function D : [— r^*, 0] — >■ M such that 

sup \u{t,x,y)-A-x^D{t)-^^y^-''\<r^+°' (8.2) 

[-r2s,0]xB+ 1-a 

Moreover, \A\ is bounded by a universal constant C depending only on dimension, s, and \ \b\\c^-2s 
and so are \D{t)\ and \D'{t)\. 

The idea of the proof is to use Lemma 17.11 from the previous section to show that u differs 
very little from a solution v to the drift-less problem with zero right hand side. This function v is 
smooth, so u is close to a smooth function and the result follows. 

Proof of Lemma [8A[ From Proposition l6.1[ the function u is Holder continuous in [—1/2, 0] x i?°/2 
and on {t} x for all t G [-1/2,0]. 

Let v be the solution to the fractional heat equation in [—1/2,0] x B^^^ with u as boundary 
values: 

vt{t,x,0) - lim y^dyvit^x^y) =0 on (-1/2,0] x B°/2 

divy^Vw in [-1/2,0] x B+2 (8.3) 

v{t,x,y) = u{t,x,y) on ([-1/2,0] x dB+) U ({-1/2} x B+^). 

From Lemma mi |u — t;| < e in [—1/4, 0] x B^^^, where e is arbitrarily small depending on the 
choice of 5. 

From Proposition 14.11 the function v is C°° respect to x, and dtv and dtVv are bounded in 
the interior of [—1/2, 0] x -Bjyj with estimates depending on the oscillation of v, which is less than 
one. In particular, \i A = Vw(0, 0, 0) and D = v{t, 0, 0), v satisfies 

\v{t,x,0)-D{t)~A-x\ < C{\x\\t\ + \xf) < Cr^+^' in [-r^^O] x B+ . 
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Also from Proposition 14. 11 D'{t) is bounded depending only on dimension and s. 
Since v satisfies the equation (|8.3p . then D'{t) ~ dtv{t, 0,0) = \imy^oy'^dyv{t,0,y). From 
Lemma [2.21 v has an expansion of the following form in [— r^^jO] x B^. 

yl-a 

v{t, X, y) = v{t, X, 0) + I lim y^dyvH, x, y) + 0(y2) 

I — a y-*-o 

= A-x + D(t) + ^^y'-'' + 0(r^+^') 
1 — a 

Meaning that the last error term 0(r^+-^*) is bounded by Cr^^'^'^ = Cr^^"' for some constant C 
independent of u and v. 

Therefore, if a < 2s, we can choose r small such that 

|^;(t, X, y) - D{t) -A-x- ^y^-^ < Cr^+^' 

1 — a 

< iri+" 
- 2 

Now we choose 5 so that the e that we obtained from Lemma [7.11 is less than ir^+". So we 
obtain 

\uit,x,y) - D(t) -A-x- ^^^y^-'^l < lr^+" + 2e < r^^" 

1 — a 2 

which finishes the proof. □ 



9 The proof of the main theorem 



We start the proof of our main theorem by proving the result in the case b{t, 0) = and f{t, 0) = 0. 
The proof of the general case follows below and consist of combining this lemma with a change of 
variables following the flow of the vector field. 



Lemma 9.1. Let a G (0,2s), and u be a solution of 

ut + b ■ Vu - lim y°-dyu{t, x, y) = f on (-1, 0] x i?° 

'y->0 

divy"Vw = m[-l,0]xB+ 
Assume u(0, 0, 0) = 0, ||wI|^ooq_i o]xb+) — ^ '^'^^ 



(9.1) 



sup 



(t,2;)G[-l,0]xSi V F 



\bit,x) 



\l-2s+a 



< S and sup 

(t,a;)e[-l,0]xi 



\f{t,x)\ 

T.|l-2s+Q 



< s 



(9.2) 



(in particular b(t, 0) = and f{t, 0) = for all t). Then there exist A G M" and D : [—1/2,0] 
such that 



u{t,x,y)-A-x-D{t)-^y'-- 

1 — a 



< C(|a-| +y + t^) 



l+a 



(9.3) 



Proof. Let p > be the radius r of Lemma \SA\ We will prove iteratively that there exists Ak and 
Dk uniformly bounded such that 



sup 

(*,a:,i/)G[-p2'>fc,0]xB\ 



u{t, x,y)-Ak-x- Dk{t) - ^^^y^-" 

1 — a 



(9.4) 
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For fc = 0, we choose A — and D = and (|9.4p holds. Now, let us assume we have shown 
(|9.4p up to some value of k and prove it for fc + 1 . 
Let 

vit, X, y) = (^uip^^H, /x, p'^y) ^p'^A.-x- Dip'^H) - p'^^'^P^^y^'^'^ 

and b = p^'^''^'^^^b{p'^'''^t, p^x) and / = p^^"^"^^^""^ f{p'^'''^t, p^x). We have that v solves 
z;f(<,a;,0) + 6- Vw- lim y''ayw(f, x, = f - p-''°'b ■ Ak on (-1,0] x B?, 

divy''Vw = in[-l,0]xB+. 



(9.5) 



Note that f — p^^°'b- Ak is small because of the assumption (|9.2p and the fact that Ak is uniformly 
boimded (independently of 5). 

Note that the term —D{t) adds a —D'{t) to the derivative of w respect to time, which is 
compensated by the term — -jz^y^"" which does the same thing with the boundary derivative 
limj^^o y°'dyv{t,x,y). 

From the inductive hypothesis, we know that < 1 in [—1,0] x Bf and from (|9.2p . we have 
that |5| < 5 and |/| < 5 in [—1,0] x . We can apply Lemma [8.11 to v to obtain that there are 
A and D uniformly bounded such that 



v{t, X, y) - t.(0, 0, 0) - A • X - D{t) - ^^y^-'^ 

\ — a 



Therefore (EH) holds for fc + 1 with Ak+i = Ak + p''°'A and Dk+i ^Dk+ p'''^"+"-^D. This finishes 
the proof of (|9.4p by induction in fc. 
Note that by construction 

i<k 

So and \Dk\ are bounded by a convergent geometric series, and then they stay bounded 
uniformly in fc. Moreover, if we set 

A = lim Ak 
D = lim Dk 

then, also by the geometric series bound, \Ak-A\ < Cp-'"^ and \Dk-D\ < Cp-^^'^+''\ Therefore, 
(|9.4p implies (|9.3p and we finish the proof. □ 

Remark 9.2. In the proof of Lemma [9.11 we see the importance of rewriting the equation (|l.ip 
using the extension as in (|2.2p . The slope with respect to time Dt is compensated by subtracting 
Dy^~°- /{\ — a) in the extended (artificial) variable. The terms fit together nicely and the error 
terms from the nonlocal operators are completely avoided. Indeed, in this supercritical case, those 
error terms would be very difficult to control without using the extension. 



Remark 9.3. If a = 0, Lemma WA\ still applies with A = 0, and all Ak vn the proof equal to zero 
as well. 
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We now have all the results needed to prove the main theorem. 

Proof of Theorem Let u be a bounded solution to (|l.ip in [—1,0] x R". By the usual trans- 
lation and dilation argument, it is enough to show that the estimate holds at the point (0, 0). We 
will prove that there is a vector A G M" such that 

\u{0,x) - u(0,0) ^A-x\< (^\\u\\l^ + "^WfWc^-.s^^j (9.6) 

where S is the constant from Lemma [9. II and r depends on ||6||pi+2s-c, ■ 
We first normalize the solution by considering 

u{r'^H,rx) -u{0,0) 



f 
b 



H\l^ + ^||/||ci-2.+o 
\u\\l=° + ^-5— ll/llci-2»+° 



b{r^''t, rx) 



In this way, ||ti||ioo < 1, and also 1 1/| |(^i-23+c, < 5 and 1 16| |(7i-2ii+a < (5 if r is small. 

Note that u, b and / also satisfy the equation (jl.ip . We would finish the proof if we could apply 
Lemma WA\ to u, b and /. We extend u and / to the upper half space. Note that the extension of 
u has the same norm as u by the maximum principle (indeed, all that is needed is that the 
extension of u is bounded in [—1,0] x Bf). The only hypothesis we are missing in order to apply 
Lemma [93] are (|9.2p . We will make another change of variables. 

We solve the following ODE backwards in time: 

1/(0) = 0, 

V{t) = h{t,V{t)) for t e (-1,0). 

Since b is continuous, the ODE has at least one solution. Moreover, |y| < C since b is bounded. 
We also define ^ 

S{t)= [ /(s,0)ds. 



which is clearly the solution of 



Therefore wc now consider 



S'(O) = 
S{t)^-f{t,()) 



u*{t,x,y) = u{t,x + V{t),y)~S{t) 
b*{t,x) = b{t,x + V{t)) - b{t,V{t)) 
f*{t,x) = fit,x)-fit,0). 

And they solve the equation 

< + 6* • Vu* - lim y^dyu* = f* in (-1, 0] x 
div y^Wu* =0 in [-1, 0] x 
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Moreover, since ||/||ci-2s+<:, < S and ||6||(;;;i-2s+c, < 6, then /* and b* satisfy (|9.2p . So, we can 
apply Lemma [93] to u* to obtain that there is A* G R" and Z?* : [— 1/2, 0] — > M such that 



"(t, X, y) - u*(0, 0, 0) - • X - Djt) - ^^*^y^-^ 

1 — a ' 



But note that u(fi, x, y) = u* (0, x, y) and 



m(0,x,j/) 

Therefore, if we set A ^ A* 



\f\\ci-2s+cj v*{0,r ^x,r' 
, we obtain 



which finishes the proof (recall that r depends on ||6||pi-2a+a). 



□ 



Proof of Theorem ] 1. SI Repeating the proof of Theorem 1 1.1 1 but with a = and A ^ 0, Theorem 
11.21 follows. The smallness of h follows from (|1.2p for small r. □ 



Remark 9.4. If we follow the iteration in the proof of Theorem 11.11 in the vanishing viscosity 
framework; we see that the artificial viscosity term eAu gets magnified in each successive scaling 
by p*^^*"^-*. For a large value of k, this value would be above the threshold to apply Lemma T?. II 
But at that point the second order diifusion takes over, and the rest of the iteration is trivial using 
the regularization of the perturbed second order heat equation. 
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